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The 3 He(e, e'p)d reaction is studied using an accurate three-nucleon bound state wave function, a 
model for the electromagnetic current operator including one- and two-body terms, and the Glauber 
approximation for the treatment of final state interactions. In contrast to earlier studies, the profile 
operator in the Glauber expansion is derived from a nucleon-nucleon scattering amplitude, which 
retains its full spin and isospin dependence and is consistent with phase-shift analyses of two-nucleon 
scattering data. The amplitude is boosted from the center-of-mass frame, where parameterizations 
for it are available, to the frame where rescattering occurs. Exact Monte Carlo methods are used 
to evaluate the relevant matrix elements of the electromagnetic current operator. The predicted 
cross section is found to be in quantitative agreement with the experimental data for values of the 
missing momentum p m in the range (0-700) MeV/c, but underestimates the data at p m ~ 1 GeV/c 
by about a factor of two. However, the longitudinal-transverse asymmetry, measured up to p m ~ 
600 MeV/c, is well reproduced by theory. A critical comparison is carried out between the results 
obtained in the present work and those of earlier studies. 

PACS numbers: 24.10.-i,25.10.+s,25.30.Dh,25.30.Fj 



I. INTRODUCTION 

Recent experiments at JLab have yielded beautiful data for the 3 He(e, e'p)d cross section and longitudinal-transverse 
asymmetry up to missing momenta p m ~ 1.1 and 0.6 GeV/c p|, respectively, and the three-body breakup cross section 
3 He(e, e'p)pn for missing momenta and energies (E m ) in the range p m up to 0.84 GeV/c and E m up to 140 MeV @|. 
These data have spurred renewed interest in these reactions, which has led to a series of papers [H, Hi El > dealing 
with the theoretical description of the proton-knockout mechanism and, in particular, with the treatment of final 
state interactions (FSI) at GeV energies. 

In the present work, we contribute to this effort. As in Refs. we adopt the Glauber approximation to describe 
the rescattering processes between the struck proton and the nucleons in the recoiling deuteron. However, in contrast 
to the studies of Refs. |E S> IE El> we re tam the full spin and isospin dependence of the nucleon-nucleon (NN) 
scattering amplitude from which the Glauber profile operator is derived, and do not make use of the factorization 
approximation, which allows one to write the (e, e'p) cross section on a nucleus in terms of a proton cross section times 
a distorted spectral function. A number of issues pertaining to the parameterization of the TV TV amplitude, and its 
boosting from the center-of-mass (cm.) frame to the frame in which rescattering occurs, are discussed in considerable 
detail. 

The spin dependence of the NN amplitudes will turn out to play an important role in the high p m region where 
double rescattering effects become dominant. It is not clear that one is justified in ignoring it, particularly in the 
analysis of experiments, such as the 4 He(e, e'p) 3 H reaction Q, in which the polarizations of the ejected proton are 
measured |8|. 

The present paper is organized as follows. In Sec. [H] we briefly review the bound-state wave functions and the 
model for the electromagnetic current operator, while in Sec. IHII the Glauber approximation, the profile operator, 
and the parameterization for the NN scattering amplitude are discussed. Next, in Sec. II VI the relevant formulae for 
the (e, e'p) cross section and their limits in the plane-wave-impulse-approximation (PWIA) are summarized, while 
in Sec. El the Monte Carlo method, as implemented in the present calculations, is described. Finally, in Sec. IVII we 
present a detailed discussion of the results, including a critical comparison between the present and earlier studies, 
while in Sec. IVIll we summarize our conclusions. 
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II. WAVE FUNCTIONS AND CURRENTS 



In this section we briefly describe the 3 He wave function and the model for the nuclear electromagnetic current. 
The discussion is rather cursory, since both these aspects of the calculations presented in this work have already been 
reviewed in considerable detail in a number of earlier papers. References to these are included below. 



A. Bound-state wave functions 



The ground states of the A— 3 nuclei are represented by variational wave functions, derived from a realistic Hamilto- 
nian consisting of the Argonne Via two-nucleon (NN) Q and Urbana-IX three-nucleon (NN N) [l Ol interactions — the 
AV18/UIX Hamiltonian model — with the correlated hyperspherical-harmonics (CHH) method The high accuracy 
of the CHH wave functions is well documented as is the quality of the AV18/UIX Hamiltonian in successfully 
and quantitatively accounting for a wide variety of three-nucleon bound-state properties and reactions, ranging from 
binding energies, charge radii, and elastic form factors [l2l fl3L IT^j to low-energy radiative and weak capture cross 
sections and polarization observables [lfij . to the quasi-elastic response in inclusive (e, e') scattering at intermediate 
energies [lfij . 



B. Electromagnetic current operator 

The nuclear electromagnetic current includes one- and two-body terms, 



p(q) = ^2pi(<l)+ ^'( q ) 

i— 1 i<j—l 
A A 

j(q) = ^^ji(q) + -Mq) 

i<j=i 



(2.1) 
(2.2) 



i=l 



The one-body current and charge operators have the form recently derived by Jeschonnek and Donnelly [l7| from an 
expansion of the covariant single- nucleon current, in which only terms dependent quadratically on the initial nucleon 
momentum (and higher order terms) are neglected. In momentum space, they are explicitly given by 



PM) = ^G E + [g m \ Ge j a, • , q s p, , 



(2-3) 



j»(q) 



TO V 



G t 



-G 



M 



2m 



-Gn 



m q x <T i 



2m 



q • (Ji q x p 4 



(2.4) 



where q (u>) is the virtual photon three-momentum (energy) transfer, Q is the four-momentum transfer with Q 2 =q 2 —uj 2 
> 0, r) is defined as rj = Q 2 /(4m 2 ), to being the nucleon mass, and Pi and <Tj are the momentum and spin operators 
of nucleon i, respectively. The nucleon Sachs form factors Ge and Gm ar e defined as 



Ge/M — ^ [(Ge p /M p + G En /Mn) + T i,z (G Ep j 'Up ~ Geu/Mu )} > (2.5) 

where the Q 2 dependence is understood, and r z is the z-component of the isospin. The Hohler parameterization 0| 
is used for the proton and neutron electric and magnetic form factors GEp/ E n an d Gmp/Mti- 

The form adopted above for the one-body currents is well suited for dealing with processes in which the energy 
transfer may be large (i.e., the ratio of four- to three-momentum transfer (Q/q) 2 is not close to one) and the initial 
momentum of the nucleon absorbing the virtual photon is small. Thus, its use is certainly justified in quasi-elastic 
kinematics for moderate values of the missing momentum. Note that in the limit (Q/q) 2 ~ 1 one recovers the standard 
non-rclativistic expressions for the impulse-approximation currents (including the spin-orbit correction to the charge 
operator). 
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The two-body charge, py(q), and current, jiy(q), operators consist of a "model-independent" part, that is con- 
structed from the TV TV interaction (the AV18 in the present case), and a "model-dependent" one, associated with the 
excitation of nucleons to A resonances (for jy only) and pir^f and mn^f transitions (for a review, see Ref. 0] and 
references therein). Improvements in the construction of the model-independent two-body currents originating from 
the momentum-dependent terms of the NN interaction have been recently reported in Ref. |l9j| . In this latter work, 
three-body currents associated with NNN interactions have also been derived. Both these refinements, however, are 
expected to have little impact on the results of the present study, and therefore are not considered any further. 

The present model for two-body charge and current operators is quite realistic at small momentum transfers. 
However, for processes involving momentum and energy transfers of order 1 GeV, such as the 3 He(e,e'p)rf reaction 
under consideration here for which q=1.50 GeV/c and w=0.84 GeV, it is likely to have additional corrections. 

III. FINAL STATE INTERACTIONS: GLAUBER APPROXIMATION 

In the kinematics considered in Sec. I VII the proton lab kinetic energies are typically of the order 0.5 GeV or larger. 
These energies are obviously beyond the range of applicability of NN interaction models, such as the AV18, which 
are constrained to reproduce NN elastic scattering data up to the pion production threshold. At higher energies, 
NN scattering becomes strongly absorptive with the opening of particle production channels. Indeed, the pp inelastic 
cross section at 0.5 GeV increases abruptly from about 2 mb to 30 mb, and remains essentially constant for energies 
up to several hundred GeV [20|. 

On the other hand, the small momentum transfer which characterizes scattering processes at high energies makes 
the Glauber approximation |2j particularly well suited in this regime. Another advantage is its reliance on a NN 
scattering amplitude, which is fitted to data. It is the approach we adopt in the present work to describe the wave 
function of the final p+(A — 1) system as 

i/i(p +( A ~Vf; GLB) = -L ]T e v G(A; 1...A— 1) &~* X Mw) e*/ "' 1;/) , (3.1) 
v A v 

where x<r(j>) represents a proton in spin state a, </> CT/ (/) denotes the wave function of the (A — l)-system with spin 
projection o~f, and Ri....a_i is the center-of-mass position vector of the A — 1 nucleons in this cluster. The sum over 
permutations V of parity e-p ensures the overall antisymmetry of ip(p GLB). 

The operator G(A; 1 . . . A— 1) inducing final-state-interactions (FSI) can be derived from an analysis of the multiple 
scattering series by requiring that the struck (fast) nucleon (nucleon A) moves in a straight-line trajectory (that is, it 
is undeflected by rescattering processes), and that the nucleons in the residual system (nucleons 1, . . . ,A — 1) act as 
fixed scattering centers [22, (the so-called frozen approximation) . It is expanded as 

A-l 

G = 1 + ^2(-) n G (n) , (3.2) 

n=l 

where G^ represents the n th rescattering term, and therefore for an A-body system up to A — 1 rescattering terms 
are generally present. The leading single-rescattering term reads 

A-l 

G ( - 1 \A;l...A-l) = Y / 8(^A)r iA (h lA ;s lA ) , (3.3) 

where Zi A and hi A denote the longitudinal and transverse components of — r A relative to p, the direction of the 
nucleon momentum, 

z iA = p • (r; - r^) , r.; -r A = b iA + z iA p , (3.4) 

and the step-function 9(x), 9{x) = 1 if x > 0, prevents the occurrence of backward scattering for the struck nucleon. 
The "profile operator" Ti A , to be discussed below, is related to the Fourier transform of the NN scattering amplitude 
at the invariant energy ^JSi A . 

The double-rescattering term, relevant for the present study of the 3 He(e, e'p)d reaction, is given by 



4 



A-l 

G^(A;1...A- 1) = S{z l] )e{z ]A )T lA {h lA -s lA )T lA {h ]A -s ]A ) , (3.5) 

where the product of 0-functions ensures the correct sequence of rescattering processes in the forward hemisphere. 
For example, if Zj — z A > and Zi — Zj > 0, then nucleon A scatters first from nucleon j and then from nucleon i. 
Note that the operators Ti A and Tj A do not generally commute. 



A. The profile operator Fij 

At this stage it is useful to specify the kinematics of the various rescattering processes occurring in the Glauber 
expansion. Specializing to the 3 He(e, e'p)d reaction of interest here, the single- and double-rescattering terms are 
illustrated schematically in Fig. ^ ln this figure, nucleon 3 denotes the knocked-out nucleon with momentum P3=p 
and energy E 3 =E in the lab frame, while nucleons 1 and 2, making up the deuteron, have momenta p! and p 2 , 
respectively, with Pi+P2=Pc; (again, in the lab frame). The black solid circle represents the AW scattering amplitude. 
In the single-rescattering case, the NN amplitudes in the two terms (panel a) of Fig. only shows one of them) are 
evaluated at the invariant energies y/sa, i=l,2, with 

s l3 = (Ei + E 3 ) 2 - ( Pi + p 3 ) 2 

~ 2m 2 + E^Jp 2 +4m 2 p p d , (3.6) 

where in the second line it has been assumed that i) the nucleons are on their mass shells, and ii) nucleons 1 and 2 
in the recoiling deuteron share its momentum equally, Pi ~ Pd/2. The momenta of nucleon 3 and nucleon i, i=l,2, 
after rescattering are p — k and Pd/2 + k, where k denotes the momentum transfer. The spectator nucleon (j ^ i) 
has momentum pd/2. Thus, the "rescattering frame"we refer to in the following is defined as that in which nucleon 
3 (nucleon i) have initial and final momenta p and p — k (pd/2 and p^/2 + k), respectively. 

In the case of double rescattering, panel b) of Fig. ^ a similar analysis can be carried out. In particular, it leads 
in the eikonal limit to the approximation in which both NN amplitudes are evaluated at the invariant energies 
y/sis ~ v^23' as obtained in Eq. I|3.6(l . 

The profile operator Fy is related to the NN scattering amplitude in the rescattering frame, denoted as Fij(k; s), 
via the Fourier transform 

r« (b; a) = ^— J d 2 k e- ik b ^- (k; s) , (3.7) 

where, in the eikonal limit, the momentum transfer k is perpendicular to p. The isospin symmetry of the strong 
interactions allows one to express Fij as 

Fij — Fij^ + + Fij-Ti ■ Tj , (3.8) 
where the are related to the physical amplitudes for pp and pn scattering via 



F ijt± = « 2 ij . (3.9) 

Available parameterizations of the pp and pn amplitudes are given in the cm. frame |22|. and therefore one needs 
to boost these from the cm. to the rescattering frame. In the discussion that follows, the general form of the 
cm. amplitude and its parameterization are described first, while the boosting procedure adopted in the present work 
is illustrated next. 
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B. The NN scattering amplitude in the cm. frame 

It is well known that the most general form for the NN scattering amplitude in the cm. frame reads (see Refs. [20LI2S 
and references therein) 



(apr^rCM) = E K N {s,tf)0™ , (3.10) 

m— 1 

where p and p' denote the initial and final nucleon momenta, respectively, the 's are functions of the invariant 

energy y/s and momentum transfer k. 2 (with k = p — p'), and a possible choice for the five operators O™, especially 
convenient for our present purposes, is that given in Ref. |22| . 

O™ = 1 , <Tj • CTj , i (<Xj + <Tj) • k x e , er.; ■ k <Xj ■ k , <j; L ■ ecrj ■ e . (3-H) 

Here the unit vector e is defined as e = (p + p')/ I P + F' I (— F/ I P I m the eikonal limit), and the overline on the 
various quantities in the equations above is to indicate that they are given in the cm. frame. The factor (2ip) _1 in 
the l.h.s. of Eq. I|3.10|l is conventional. 

— NN 

The functions F m are parameterized as 



F;"( 8| k-)=C(«)exp -/C(*)k 2 l , (3.12) 



-NN 



where the (s) and (3 m (s) coefficients depend on s and are generally complex; in particular, the forward, spin- 
independent amplitude F 1 (s, 0) is given by F 1 (s, 0) = a (1 — ip NN ) / (87r) (an invariant quantity), where a NN 

is the total cross section and p NN is the ratio of the ima gina ry to the real part of F 1 (s, 0). In the present work 

these coefficients are taken from Tables III and IV of Ref. 22]: they were obtained by fitting NN amplitudes derived 

from the phase-shift analysis of the VPI group |24|. The tabulations are for s values in the range (3.92-5.09) GeV, 

corresponding to lab kinetic energies (210-831) MeV, for both pp and pn scattering. A detailed assessment of the 

accuracy and limitations of these parameterizations (of course, in relation to the NN scattering data available up 

to 1981) can be found in Ref. j^. Here, it suffices to note only that they are reasonably accurate for the central 

2 

and single spin-flip terms for k up to ~ 0.1 (GeV/c) 2 ; however, double spin-flip terms are not well fitted by (single) 
Gaussian functions. 

It would be desirable to update and improve the parameterizations of Ref. p2| by using phase-shift analyses based 
on the current NN database. Unfortunately, at the higher lab energies of interest here (say, above 500 MeV) the 
increase in the size of this database has been rather modest since the early 1980's. 

It is worth pointing out that all previous Glauber calculations of A(e, e'p) processes we are aware of have ignored 
the double spin-flip terms in the scattering amplitude, corresponding to m =2, 4, and 5 in Eq. I|3.10|l : indeed, most 
have only included the spin independent term (see, for example, Refs. |3l H. l23l l2q ). 

Finally, note that, if the amplitude in Eq. I|3.10l) were to consist only of the scalar term (m=l), then the trans- 
formation to the lab frame would be unnecessary, since (2ip) _1 F ^ would be an invariant function of s and the 

2 

four-momentum transfer squared t=— k . However, the presence of the spin-dependent terms spoils this simplicity. 



C. Boosting the NN scattering amplitude to the rescattering frame 

In order to boost the NN scattering amplitude from the cm. to the rescattering frame, we adopt the procedure 
described in Ref. |26| , although its practical implementation in the present work is approximate for reasons discussed 
below. First, an invariant representation of the amplitude is introduced in terms of scalar, vector (7^), tensor (c^), 
pseudoscalar (7 s ), and axial vector (7 5 7 M ) combinations of Dirac matrices: 



NN 



E 



(3.13) 
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where the five operators A™ are defined as 

Am- 1,...,5 -i u ai> 5 5 5 /^5 /o i a\ 

ij = 1 > 7f 7j> , < o-j,/*!/ . 7* 7j , 7i 7i 7i7j> ■ (3-14) 

The relation between the functions J 7 ^ 1 ^ and F m in the cm. frame follows by noting that 

^(PO^(-PO^^(P)M-P) = xlxl. [(^Pr^^iKs)] x„Xa 3 , (3.15) 

where the u a are (positive-energy) Dirac spinors with u a = u^j , and x<x are two-component Pauli spinors. This 
leads in the eikonal limit p ~ p ' to 

5 

T% N (s,k 2 ) = ^M mn (p,k 2 )Ff N ( S ,k 2 ) , (3.16) 

n=l 

where the matrix M is the inverse of that given in Tables I and II of Ref. [2^|. In fact, the k dependence in the 
matrix M is neglected, i.e. M = M (p, 0), in the present work, which allows one to write , m = 1, . . . , 5, as a linear 
combination of terms, each having, as far as the momentum transfer dependence is concerned, the same Gaussian 

functional form as F n . This turns out to be convenient when performing the Fourier transform in Eq. 

Next, having determined the functions > the scattering amplitude in the rescattering frame is obtained from 

xlxl, [(21P)- 1 F» N (k, s)] XatXvj = u a , (p - k)u a , (p d /2 + k)T^ N u ai (p) U(7j (p d /2) , (3.17) 

where in practice the dependence upon p^/2 in the spinors of particle j has been neglected (in this limit, the rescat- 
tering and lab frames coincide). This approximation is justified at low pd (corresponding to low missing momenta), 
but is clearly unsatisfactory at high pd (pd — 1 GeV/c in the experiment of Ref. I}). The resulting F^ N (k, s) has 
central, single and double spin-flip terms. 

The approximations above — neglecting the dependence of the matrix M on the momentum transfer in Eq. I|3.16|l . 
and the momentum p^/2 in Eq. (|3.17|l -made in boosting the NN amplitude from the cm. to the rescattering frame, 
have been dictated by computational convenience rather than by necessity, and could be removed. This latter task, 
however, is beyond the scope of the present work. 

IV. CROSS SECTION AND RESPONSE FUNCTIONS 

To set the stage for the discussions that follow in later sections, it is helpful to give the expression of the five-fold 
differential cross section for the A i(e, e'p)*-" 4-1 ^/ process (for a derivation, see Ref. |27|'): 

d 5 

, rwjo, no = P ^ CTMott /rec-p "FT" [vlRl + V T Rt + V LT RltCOs{4>) + VfT RtT COS (20)] , (4.1) 

aF' e d\V e a\l ki hf 

where E' e is the energy of the final electron, Vl' e and VL are the solid angles of, respectively, the final electron and 
knocked-out proton, nif is the rest mass of the (^4-l)-cluster (assumed bound here), p and E (p/ and Ej) are the 
momentum and energy of the proton ((A-l)-cluster) in the lab system, <j> is the angle between the electron scattering 
plane and the plane defined by the momenta q and p, and the recoil factor / rec , or rather its inverse, is defined as 



r 1 = 

J rec 



(4.2) 



The Mott cross section UMott and the coefficients v a , a~L, T, LT, and TT, are defined in terms of the electron 
kinematical variables, as given in Eqs. (2.19) and (2.27a)-(2.27rf) of Ref. (note, however, that in that work Q 2 is 
taken to be negative). 
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The response functions R a involve matrix elements of the charge and current operators between the initial i and 
final p+^ A ^^f nuclear states, and depend on the momenta q and p, the angle 9 between them, and the energy transfer 
lu. In a schematic notation, they are given by 

Rl = \{p+^f\p(qz) \ A l)\ 2 , (4.3) 

Rt = \(p+ {A - 1] f\lx{qi)\ A i)\ 2 , (4.4) 

R LT = 2V2 {p+ {A -^f | P {qi) \ A i)* {p+ {A - 1] f | j x (qz) \ A i) , (4.5) 

Rtt = -\{p+^f\U<li)\ A ^)\ 2 + \{P+ {A - 1) f\jy{<li)\ A f}\ 2 , (4.6) 

where the z-axis has been taken along q, which also defines the spin-quantization axis, jj_ denotes the components of 
the current transverse to q, and the average over the initial, and sum over the final, spin projections are understood. 

A. The plane-wave-impulse-approximation 

In the plane- wave-impulse-approximation (PWIA) limit, in which final-state interactions (FSI) effects between the 
knocked-out proton and the nucleons in the recoiling cluster are ignored (i.e., the operator G in Eq. (|3.1|1 is set to one), 
the response functions i?^ WIA can be expressed, neglecting two-body terms in the electromagnetic current operator, 
as 

Rl WIA = r%N pf (p m ), (4.7) 

where the r v a denote appropriate combinations of kinematical factors with the proton electric and magnetic form 
factors — those corresponding to the one-body currents of Sec. Ill Bl are listed in Appendix lAl — and N p f(p m ) is the 
p+(A-l)-cluster momentum distribution, defined as 



N, 



pf 



{Pm) = 2j~Ti £ KWp-)I 2 ' ( 4 - 8 ) 



with 



A p J,<r f -, at (Pm) = ^ 7^)3 e -W<rA-Ri...A-0 % t {A;p) ^ M....1 ! : / : , -, T ; 1 . . . ,1: / 1 . (4.9) 

Here Jj is the total angular momentum of the initial state, and p m is the so-called missing momentum, p m =— P/=p q. 
The normalization integral 



N f = 4tt / d Pm pi N pf (p m ) (4.10) 
Jo 

gives the number of / clusters in the ground state A i [2^|. In 3 He the number of deuterons is calculated to be 
about 1.34 (see below), which implies that in 3 He (to the extent that it is a pure total isospin T=l/2 state), out 
of a possible number of 1.5 pairs of nucleons in isospin T=0 states, roughly 90% of them are in the deuteron state. 
Similarly, in 4 He the number of tritons is found to be ~ 1.68, and so in 4 He (again, ignoring admixtures of states with 
T > 0, induced by small isospin symmetry-breaking interactions), about 85% of the nnp clusters are in the triton 
state [iliill^. 

The Npd(p m ) momentum distribution, obtained with CHH wave functions corresponding to the AV18/UIX Hamil- 
tonian model, is shown in Fig. up to missing momenta of ~ 1 GeV/c. A number of realistic interactions are 
currently available, such as, for example, the CD-Bonn [3l| or Nijmegen y^] NN and Tucson-Melbourne NNN 
interactions, and therefore the question arises of how sensitive to the input Hamiltonian are the high momentum 
components of this momentum distribution. This issue is especially relevant here, since the kinematics of the JLab 
(e, e'p) experiments cover a broad range of p m values, as high as 1.1 GeV/c. It is addressed in Fig. |3 where pd 
momentum distributions, obtained with various combinations of two- and three-nucleon interactions, are compared 
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with each other up to p m ~ 1 GeV/c. All results, but for those labeled AV18/UIX-CHH, are obtained with Faddeev 
wave functions [34j . 

A couple of comments are now in order. First, in the p m range (400-800) MeV/c, there is a significant model 
dependence: the CD-Bonn N p d(p m ) is about a factor of 2 smaller than the AV18 one. This is likely a consequence of 
the fact that the tensor force is weaker in the CD-Bonn than in the AV18. The pd overlap in 3 He has S- and D-state 
components, and the associated D-state contribution to N p d(pm) indeed becomes dominant at p m ~ 400 MeV/c, it is 
responsible for the change of slope in N p d(p m ). 

Second, the Faddeev and CHH wave functions corresponding to the AV18/UIX Hamiltonian model lead to pd 
momentum distributions, that are slightly different only in the p m region around 400 MeV/c, where the S-wave 
contribution changes sign. 



V. CALCULATION 



Nuclear wave functions, for an assigned spatial configuration R = (rj, . . . , r^), are expanded on a basis of 
K=2 A A\/[Z\(A — Z)\] spin-isospin states for A nucleons (Z is the number of protons) as 



K 

^(R)=5> fc (R)|*c), (5.1) 

fc=i 

where the components ^fc(R) are generally complex functions of R, and, in the case A=3 and Z=2 as an example, 
the basis states \ k)= \ (n l)i(p D2& Ih) , I (nt)i(p|)2(pj.)3), ■•• , | (pt)i, (pih{n ■■■ ■ Matrix elements of the 
electromagnetic current operator are written schematically as 




where [Ojy(R)] denotes the matrix representing in configuration space any of the one- or two-body charge/current 
operators. Matrix multiplications in the spin-isospin space are performed exactly with the techniques developed in 
Ref. (3^. The problem is reduced to the evaluation of the spatial integrals, which is efficiently carried out with Monte 
Carlo (MC) methods, although these are implemented differently in the present study than they have been in the 
past |28ll3q |. 

To illustrate these methods, consider the PWIA calculation of the one-body charge operator in the process 
3 Ee(e,e'p)d: 



P PWIA ~]T / dRe- i P- r3 e- i P/- Rl2 Vfe(ri-r 2 ;/)e i< J- r3 p fe i(q;-iV 3 )Vi(R;0 , (5.3) 

k,l J 

where ip(i) and ip(f) denote the 3 He and p+d cluster wave functions, respectively, and Pfcz(q;— 1V3) is kl matrix 
element of the charge operator in Eq. I|2.3|) . Its dependence on the momentum operator, through the spin-orbit term, 
has been made explicit. It is convenient to introduce the Jacobi variables, x=r 3 — R12, y=i"i — y%, so that the integral 
reads 



p pwia dxdye- i P- x ^(y;/)p fc ; (q;-4 Vx ) lMx,y;») , (5.4) 

where p m is the missing momentum. In Refs. 0[3(J, this integral was evaluated with the Metropolis algorithm |37j 
by sampling the coordinates (x, y) from a probability density, taken to consist of products of central correlation 
functions from VMC wave functions (see Ref. |2^). While this procedure is satisfactory at low p m , in the sense 
that statistical errors are small, it becomes impractical as p m increases, due to the rapidly oscillating nature of the 
integrand. Indeed, this problem is evident in the VMC calculation of the pd and pt momentum distributions, shown 
in Ref. [2^] : even though the random walk consists of a large number (of the order of several hundreds of thousands) 
of configurations, oscillations in the central values persist at large p m (> 600 MeV/c). 
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In the present work, we carry out the multidimensional integrations by a combination of MC and standard quadra- 
tures, namely we write 



where the c's denote configurations (0x>y) (total number N c ), sampled with the Metropolis algorithm from a proba- 
bility density W (normalized to one), given by 



W*,y) = ^£|^(y;/)| 2 • (5.6) 

Note that W is uniform in < cj) x < 2tt. For each configuration c, the function F is obtained by Gaussian integrations 
over the x and cos^ variables, i.e. 



F(c) = E J ^ose x ) dxx 2 e'^r k {y;f)Pki (q; -^ Vx ) ^(x,y;i) . (5.7) 

As a result, the statistical errors are very significantly reduced. In Fig. [21 the CHH calculation of N p d(p m ) is carried 
out with this method, it uses a random walk consisting only of 20,000 configurations. However, convergence in the 
(x^cosQx) integrations requires of the order of (70,50) Gaussian points at the highest p m , and so the present method 
turns out to be computationally more time-consuming than the earlier version at high p m . 

Additional refinements in the present MC implementation are i) the application of gradient operators on the left, 
rather than right, wave function, and ii) the use of block averaging for a more realistic estimation of the statistical 
errors. 

Gradient operators, such as those occurring in the one-body electromagnetic current, are discretized as 



. <K...r i + 6e a ...)-1,{...r i -6e a ...) 
I d 

where 5 is a small increment ((5=0.0005 fm in the calculations reported here) and e a is a unit vector in the a-direction. 
Therefore, again in the context of the PWIA calculation above, the V x , when operating on the left, only acts on the 
plane wave in fact, an eigenfunction of V x • This further reduces statistical errors, and also ensures that the PWIA 
relations in Eq. (|4.7|l are satisfied — modulo tiny discretization errors of order (5p m ) 2 — at each configuration in the 
random walk, which would not be the case if the gradient were left to operate on the 3 He wave function to the right. 
Of course, the eigenfunction property above is spoiled, when final state interactions are taken into account; however, 
the error-reduction benefits remain. The disadvantage of the procedure just outlined is that it leads to an increase in 
computational time, since the various gradients have to be evaluated, rather than once (when acting to the right), as 
many times as the number of kinematics being considered in the calculation. 
A crude estimate of the MC error is obtained as 



A(F) 



1 iVc r 1 n c 

— YF 2 (c)~ —Yf(c) 



1/2 



(5.9) 



it assumes that the distribution F(c) is Gaussian, whereas in practice this is generally not the case. A better estimate, 
adopted in the present work, is obtained by dividing the set of N c samples into M c blocks containing n c samples each: 



^ m n c 

fm = / 



F(c) 



c=(m—l) n c + l 



(5.10) 



Then, 
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A(/) 



1 Ma I ^ Mc 

rn — 1 \ rn — 1 



1/2 



(5.11) 



where n c can be chosen large enough so as to make the distribution of f m Gaussian (in practice, n c has been taken 
of the order of 100). 



VI. RESULTS 

The predicted cross sections are compared with the data taken at JLab (E89-044) Q in Figs.Oland0] The in-plane 
measurements were carried out in quasi-elastic kinematics with the proton being detected on cither side of the three- 
momentum transfer q (the kinematics in Figs. and 01 have, respectively, 0=180 deg and deg, in the notation of 
Sec. IIVII . From these cross sections, the longitudinal-transverse asymmetry Alt is obtained as 

. cr(0 = deg) - a((f> = 180 deg) 

Alt — — 



Odeg) + cr(0 = 18Odeg) 
vlt Rlt 



vl Rl + vt Rt + vtt Rtt 



(6.1) 



where cr(0) stands for the five-fold differential cross section in Eq. I|4.1|l . Thus, the asymmetry, shown in Fig. [S] is 
proportional to the Rlt response function. Note that vlt is negative, as defined in Ref. |27j . 

In the figures, the results of three different calculations are displayed. The curve labeled PWIA is a plane-wave- 
impulse-approximation calculation, based on the CHH 3 He wave function corresponding to the AV18/UIX Hamiltonian 
model (the resulting pd momentum distribution is shown in Fig. |2J) . It neglects final-state-interaction (FSI) effects 
and contributions from two-body currents (MEC). The present PWIA results are in agreement with those of recent 
studies 0, 0: they overpredict the measured cross section for values of the missing momentum p m up to ~ 300 
MeV/c, while underpredicting it at high p nl . This underprediction is particularly severe for the kinematics with 0=0 
deg (Fig. the contribution proportional to Rlt in Eq. 1)4. 1|) is comparable in magnitude for p m > 300 MeV/c, but 
of opposite sign, relative to that from Rl and Rt (the Rtt contribution is negligible). Of course, when 0=180 deg, 
they have the same sign. 

The curve labeled GLB (GLB+MEC) shows the results of calculations using the Glauber approximation to describe 
FSI and one-body (one- plus two-body) currents. Both single- and double- rescattering terms are retained in the 
Glauber treatment of the final p+d scattering state. The profile operator is obtained from the full NN scattering 
amplitude, boosted from the cm. to the rescattering frame as discussed in Sec. IIII Gl The invariant energy y/s at 
which the rescattering process occurs, Eq. <|3.6[l . changes little over the whole range of missing momenta, it corresponds 
to a lab kinetic energy of ~ 830 MeV. Thus the parameters used for the central, single and double spin-flip terms of 
the pp and pn amplitudes are taken from the last row of, respectively, Tables III and IV of Ref. [23]. However, the 
parameters for the B pn (k) and E pn (k) terms (in the notation of Ref. corresponding to m—2 and 5 in Eq. I|3.10|l . 

are replaced by those for B pp (k) and E pp (k), since the real parts of fig 1 and 0^ are negative. 

There is satisfactory agreement between theory and experiment up to missing momenta of 700 MeV/c. At values 
of p m ~ 1 GeV/c, however, the theoretical results are smaller than the experimental values by about a factor of 
two, although they do reproduce the flattening of the cross section, as function of p m , seen in the data. Final-state- 
interaction effects play a crucial role: they reduce the PWIA cross section at low p m (< 300 MeV/c), and increase 
it very substantially at larger values of p nl . Two-body current contributions, while not large, are not negligible, 
and improve the agreement between theory and experiment. This is most clearly illustrated in the case of the Alt 
observable, Fig. [SJ 

It would be interesting to investigate the model dependence due to the input Hamiltonian used to generate the 
bound-state wave functions. We do not expect it to be large, even for p nl = (400-800) MeV/c where the calculated 
pd momentum distributions can differ by as much as a factor of two, see Fig. The cross section in this range of 
p m values results from strength shifted by FSI from the low p m region, p m up to ~ 300 MeV/c, where the N p d(p m ) 
model dependence is negligible. Clearly, a direct calculation is needed to verify whether this expectation is justified. 

The next set of figures, Figs.lBTHOl is meant to illustrate the effect of various approximation schemes in the Glauber 
treatment of FSI. In all calculations MEC contributions are included. In Figs. and [7| the results of calculations 
using i) only the central part of the scattering amplitude (curves labeled "central F only" ) and ii) the full scattering 
amplitude but neglecting boost corrections (curves labeled "full F cm ") are compared with the baseline GLB+MEC 
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results of Figs. and H (labeled here as "full F lab ")- For values of p m > 600 MeV/c, the spin-dependence of the 
NN scattering amplitude, which is ignored in all Glauber calculations of (e, e'p) reactions off few-body nuclei wc are 
aware of (see, for example, Refs. USESE^Dj leads to a very substantial increase of the cross section obtained when 
using the central term only. On the other hand, boost corrections, which in the present work are only accounted for 
approximately (see Sec. IIII Cty . seem to be small, although not negligible. 

In Figs. |H1 and the results of calculations using the full scattering amplitude in the cm. frame but including 
only single rescattering in the Glauber expansion (curve labeled "GLB-1 full F cm "), i.e. the term G (1) of Eq. (|Ol . 
are compared with the predictions obtained by retaining both single and double rescatterings (the curve labeled 
"GLB-(l+2) full F cm "is the same as "full F cm "in the previous two figures). Also shown are the PWIA results. At 
low missing momenta (below 300 McV/c) interference between the plane- wave and single-rescattering amplitudes 
reduces the PWIA cross section. In the p m range ~ (300 — 800) MeV/c destructive interference occurs between the 
leading single- and double-rescattering amplitudes, resulting in a reduction of the cross section obtained in the GLB-1 
calculation. At the highest values of p m , double-rescattering processes are dominant. The interference pattern among 
these various amplitudes is consistent with that obtained in the calculations of Ref. Q . 

Lastly, the Alt asymmetry is found to be relatively insensitive to the various approximation schemes discussed 
above in the region where measurements are available (p m < 600 MeV/c). All calculations reproduce quite well the 
oscillating behavior of the Alt data. 

As already mentioned, most Glauber calculations of A(e, e'p) reactions have only used the central part of the NN 
scattering amplitude F^ N , i.e. 



(Hp)" 1 F.!J N (k, s) — ^ NN (s) [1 - ip NN (s)} exp [-(3 NN (s)^} , (6.2) 

where a NN (s) and p NN (s) are, respectively, the total cross section and ratio of the real to imaginary part of the 
scattering amplitude at the invariant energy yfs. While the values for a NN (s) and, to a less extent, p NN (s) are well 
known for both pp and pn over a wide range of t/s (see, for example, Ref. |25| and references therein), this is not 
the case for /3 (s). This parameter is determined by fitting either the elastic differential cross section at forward 
scattering angles (as in Ref. [2{| and references therein), i.e. da^ N (s,t)/dt for small four-momentum transfer t, or 
rather the central term of the scattering amplitude derived from phase-shift analyses (as in Ref. |22j and references 
therein). It should be emphasized that the first procedure tacitly assumes that the contributions to dcr^ N (s, t)/dt 
due to spin-dependent terms in are negligible. It is not obvious that this assumption is justified. For example, 

large cross section differences, Aa^ N , between spin orientations parallel and antiparallel to the beam direction are 
observed in pp and pn scattering |2£|. Also observed are substantial, although less dramatic, cross section differences, 
Aa^ N , between parallel and antiparallel transverse spin orientations |20j |. 

The sensitivity of the 3 He(e, e'p)d cross section, calculated using these two different values for the (3 PP parameter 
(the isospin dependence is ignored in the following discussion), is illustrated in Figs. 1111 and 1121 The curves labeled 
"central F pp : small /3 pp "was obtained with f3 pp — 0.095 fm 2 , as reported in Ref. [23 (for pp scattering). The curve 
labeled "central F pp : large /3 pp "was obtained, instead, with /3 PP =0.157 fm 2 , a value in line with that inferred from 
phase-shift analyses [22| ■ Incidentally, this large f3 pp was used in a set of unpublished calculations carried out by the 
present authors in 2004 (and referred to in Refs. 0,0,0)- 

The results of the calculation with a small /3 PP , including only one-body currents, are in quantitative agreement 
with those of Ref. Q, although at p m ~ 400 — 700 MeV/c the present predictions are somewhat larger than obtained in 
Ref. Q . These relatively small differences are presumably due to i) the breakdown of the factorization approximation 
employed in Ref. and ii) the use of the CC1 parameterization (3^ adopted in Ref. Q, rather than free nucleon 
form factors as in the present study. 

In the high p m (~ 1 GeV/c) the small f3 NN results lead to a large increase of FSI contributions. The profile function 
corresponding to a central scattering amplitude reads 



N 



(b;s) 



1 

8tt(3 nn (s) 



00 [1 



„NN 



(s)] exp 



(6.3) 



and a small /3 N , while reducing its range, makes the value at zero impact parameter larger, which leads to the FSI 
enhancement mentioned above. 

Lastly, it is interesting to note that the results of the calculation based on the full scattering amplitude derived from 
phase-shift analyses are close to those using the central term only in the amplitude with a (3 NN parameter obtained 
from the small t slope of the elastic cross section. 
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VII. CONCLUSIONS 

In this work we carried out calculations of the 3 TLe(e,e'p)d cross section for the kinematics of JLab experiment 
E89-044, spanning the missing momentum range (0-1.1) GeV/c. Final state interactions were treated in the Glauber 
approximation, including both single and double rescattering terms. In contrast to earlier studies of the same pro- 
cess 0, , the profile operator retained the full spin and isospin dependence of the underlying N N scattering ampli- 
tudes. Parameterizations for these were derived from phase-shift analyses in Ref. 22]. It would be desirable to update 
and improve these parameterizations, although the paucity of additional TV TV scattering data at energies beyond 500 
MeV collected in the last two decades or so (due in part to the termination of iViV programs at facilities such as 
Saturne and LAMPF) will presumably not alter them significantly. 

Corrections arising from boosting the amplitude from the cm. to the relevant frame for the rescattering processes 
were found to be relatively small. However, the boosting procedure was only implemented approximately. 

Theory and experiment are in quantitative agreement for missing momenta in the range (0-700) MeV/c. Rescat- 
tering effects play a crucial role over the whole range of p m , in particular double rescattering processes are responsible 
for the increase of the cross section at p m ~ 1 GeV/c, which, nevertheless, is still underpredicted by theory by about 
a factor of two. However, the flattening of the data in this p m region is well reproduced. Two-body current contribu- 
tions are relatively small, but helpful in bringing theoretical predictions for the Alt observable in significantly better 
agreement with experiment. Spin-dependent terms in the scattering amplitude are important at high p m . 

A generalized Glauber approach has been recently developed for (e, e'p) processes (sjj, El which attempts to partially 
remove the frozen approximation implicit in the original derivation |2l| . The resulting correction leads, in essence, 
to a modification of the profile operator by a phase factor — > exp(iAoZjj) T^ , with Ao fixed by the kinematics, 
Aq = LoE m /q. It was found to be numerically very small, at the most 10% at p m ~ 1 Gev/c. It may play a more 
prominent role, however, in the three-body electrodisintegration 3 He(e, e'p)pn at high missing energies. 

Glauber calculations using only the central part of the scattering amplitude with a f} NN obtained by fitting the 
low t slope of do-^ N /dt are close to those using the full amplitude derived from phase-shift analyses. As argued in 
the previous section, however, it is not clear that one is justified in ignoring the spin dependence of the amplitude 
in view of the larg e cross section differences observed in TV TV scattering involving polarized beam and target (see, for 
example, Ref. |20|). 

Future work aims at: i) extending the present Glauber approximation, based on a spin- and isospin-dependent 
scattering amplitude, to treat the p+ 3 H electrodisintegration of 4 He and, in particular, the reaction 4 Hc(e, e'p ) 3 H, 
both of which have been recently measured at JLab p],|41| (the polarization parameters measured in the 4 He(e, e' p) 3 H 
reaction have already been found to be in agreement with the results of a calculation using an optical potential |42|'): 
ii) investigating the model dependence of the present predictions for the 3 He(e, e'p)d cross section upon the input 
(non-relativistic) Hamiltonian adopted to generate the bound-state wave functions; iii) exploring the extent to which 
the use of a relativistic Hamiltonian to generate these wave functions alters the present predictions, particularly 
at high p m (the research projects in items ii) and iii) will be made possible by very recent developments of the 
hyperspherical-harmonics method in momentum space [Sj): and iv) applying the methods developed here to the 
three-body electrodisintegration of 3 He. Studies along these lines are being vigorously pursued. 
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APPENDIX A 

We list here explicit expressions for the single-proton response functions r v a defined in Eq. 14.71) : 
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where p± is (the magnitude of) the component of the proton momentum transverse to the momentum transfer q, and 
rj=Q 2 1 (Am 2 ) . The dependence of the proton electric (Ge p ) and magnetic (Gm p ) form factors on Q 2 is understood. 
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a) 



b) 



FIG. 1: Schematic illustration of single, panel a), and double, panel b), rescattering processes. Dashed (solid) lines represent 
virtual photons (nucleons), while black solid circles represent NN scattering amplitudes. 



10' 



10 L 



10" 
10" 2 
"^10" 3 

!xio- 4 

CL 

Z 10" 6 
10" 7 
10" 8 
10" 9 



AV18/UIX-CHH 

AV18/UIX 

AV18/TN99' 

CDBONN/UIX 

CDBONN/TM99' 




400 600 800 

p m (MeV/c) 



1200 



FIG. 2: The pd momentum distribution, obtained with a correlated-hyperspherical-harmonics (CHH) wave function corre- 
sponding to the AV18/UIX Hamiltonian, is compared to those obtained with Faddeev wave functions corresponding to different 
combinations of two- and three-nucleon interactions. 
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FIG. 3: The experimental data for the 3 He(e, e'p)d cross section at 0=180 deg are compared to the results of calculations 
in plane- wave-impulse-approximation (PWIA), or using the Glauber approximation without (GLB) and with (GLB+MEC) 
inclusion of two-body currents. The profile operator in the Glauber expansion is derived from the full AW scattering amplitude, 
boosted from the cm. to the the rescattering (i.e., lab) frame. 




FIG. 4: Same as in Fig. El but at <j>=0 deg. 
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FIG. 5: Same as in Fig. but for the longitudinal-transverse asymmetry rather than the cross section. 
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FIG. 6: The experimental data for the 3 He(e,e'p)d cross section at <^>=180 deg are compared to results of Glauber calculations, 
using a variety of approximation schemes for the AW scattering amplitude. See text for an explanation of the notation. 
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FIG. 7: Same as in Fig. © but at (f>=Q deg. 
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FIG. 8: The experimental data for the 3 He(e, e'p)d cross section at 0=180 deg are compared to Glauber calculations including 
only single (GLB-1) or both single and double [GLB(l+2)] rescattering terms. Boost corrections in the NN scattering amplitude 
are neglected. Note that two-body current contributions are included. 
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FIG. 9: Same as in Fig. El but at <j>=0 deg. 
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FIG. 10: The experimental data for the longitudinal-transverse asymmetry in the 3 He(e, e'p)d reaction are compared to the 
results of Glauber calculations including only single (GLB-1) or both single and double [GLB-(l+2)] rescattering terms. The 
GLB-(l+2) calculations use the full NN scattering amplitude, with (full Fi&b) or without (full F CU i) boost corrections, or its 
central term only (central F only). The GLB-1 calculation uses the full AW scattering amplitude, but ignores boost effects. 
Note that two-body current contributions are included in all calculations. 
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FIG. 11: The experimental data for the 3 He(e, e'p)d cross section at 0=180 deg are compared to results of Glauber calculations 
using the full NN scattering amplitude (curve labeled "full -Fi a t,"), or its central term only with a f3 pp parameter consistent 
with either Aa PP /At at small t (curve labeled "central F pp : small /3 PP ") or phase-shift analyses (curve labeled "central F pp : 
large /3 PP "). All curves include two-body current contributions. For the "central F pp : small f3 pp " calculation also shown are the 
results obtained with one-body currents. 
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FIG. 12: Same as in Fig. [TT] but at 4>=0 deg. 



